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I report on a calculation of the inclusive Higgs boson production cross section at hadron colliders at next-to- 
next-to-leading order in QCD. The result is computed as an expansion about the threshold region. By continuing 
the expansion to very high order, we map the result onto basis functions and obtain the result in closed analytic 
form. 



1. Introduction 

At the LHC, gluon fusion will be the most 
important mechanism for Higgs boson produc- 
tion and discovery for masses below ~ 700 GeV. 
The discovery of the Higgs boson and the subse- 
quent study of its properties therefore relies on 
a solid theoretical understanding of the gluon fu- 
sion production mechanism. Unfortunately, next- 
to-leading order (NLO) studies of inclusive 
Higgs production do not provide this solid under- 
standing. The NLO corrections are so large (of 
order 70 — 100%) that one cannot assume that 
they provide a reliable estimate of the total cross 
section. The unsettled nature of such an impor- 
tant signal clearly calls for a renewed effort to 
bring this process under control. 

Earlier this year, we computed the full NNLO 
corrections to the hadronic cross section for Higgs 
boson production Q as an expansion about the 
threshold region. The expansion was carried out 
to very high (18th!) order so that any uncertainty 
due to uncalculated higher terms would be very 
small. Very recently, this calculation has been 
confirmed by an exact NNLO calculation of the 
partonic cross section || . Shortly before this con- 
ference began, we extended our calculation of the 
expansion to sufficiently high order as to allow us 
to invert the series and obtain the exact result for 
the hadronic cross section. In addition, we have 
used the same methods to compute the NNLO 



* Supported by the U.S. Department of Energy under Con- 
tract No. DE-AC02-98CH10886. 



corrections to pseudo-scalar Higgs production, for 
which we present preliminary results Q] . 

2. The Calculation 

In the limit that all quark masses except that 
of the top quark vanish, gluons couple to Higgs 
only via top quark loops. This coupling can be 
approximated by an effective Lagrangian j^] cor- 
responding to the limit m t — > oo, which is valid 
for a large range of Mh, including the currently 
favored region between 100 and 200 GeV. The 
effective Lagrangian is 
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where G°„ is the gluon field strength tensor, H 
is the Higgs field, v « 246 GeV is the vacuum 
expectation value of the Higgs field and Ci(a s ) is 
the Wilson coefficient, which for this calculation 
we need to order (a^) j6). 

For the pseudo-scalar Higgs, we again assume 
that gluon fusion through top-quark loops is the 
dominant production mechanism. This assump- 
tion is not valid in all models, especially when the 
ratio of vacuum expectation values is large. The 
effective Lagrangian for pseudo-scalar production 

1 \di 6 X + C 2 Oal , 
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where A is the pseudo-scalar Higgs field, g t is a 
model-dependent coupling constant and G" is 
the dual of the field strength tensor: 
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The calculation breaks down into four contribu- 
tions: virtual corrections to two loops, single-real- 
emission to one loop, double-real-emission at tree- 
level and mass factorization. The virtual, single- 
real and mass factorization terms are computed 
exactly in closed analytic form. The double- 
real contribution is by far the hardest part of 
the calculation because of the complicated phase 
space integrals. We have computed this contri- 
bution by expanding the phase space integration 
about the threshold limit, where the partonic 
center-of-mass energy is close to the Higgs mass 
{Mfj/s = x — ► 1). Because Higgs production is 
dominated by threshold corrections, the series ex- 
pansion converges quite rapidly. So, we compute 
double-real-emission, and thus the partonic cross 
section, as an expansion in (1 — x) and ln(l — x). 
Note that if all coefficients are computed, this is 
an exact expression. In fact, one can obtain the 
exact result for the partonic cross section from a 
finite number of terms. 

3. Inverting the series 

If one were to know the basis functions that 
make up the exact result and one could expand 
the series out to enough terms, one could invert 
the series, mapping it onto the basis functions. 
One can obtain a reasonable ansatz for the ba- 
sis functions by examining the result for Drell- 
Yan production Using polylogagrithm iden- 
tities ^] this result can be expressed in terms of 
functions which are analytic in (1 — x) and pow- 
ers of ln(l — x) times functions which are analytic 
in (1 — x). Each of these basis functions can be 
multiplied by a pre-factor. Again using the Drell- 
Yan result as a guide and knowledge of the gluon 
splitting function, one can make an ansatz of the 
possible pre-factors. The ansatze for pre-factors 
and basis functions are shown in Table |l|. 

One sees that there are 7 pre-factors and 14 
functions meaning that if one can expand the se- 
ries result out to 98 terms, one can map the re- 
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suit onto these 98 functions. The mapping can be 
verified by computing still higher terms and com- 
paring to the expansion of the mapped functions. 
To carry out this program, we have computed the 
double-real radiation terms out to order (1— a;) 100 . 
Since the series starts at order (1 — x)^ 1 , this gives 
us 102 terms. 

It turns out that the 1/x pre-factor never oc- 
curs in the cross section (in Ref. || it appears in 
the Drell-Yan correction term, not the Drell-Yan 
cross section) and the last two basis functions, 
Li 3 ((l-x)/(l + x)) and Li 3 (-(1 -x)/(l + x)) al- 
ways occur together as the difference. Thus, after 
the fact, we see that 78 terms would have sufficed 
to determine the functional form. The additional 
24 terms significantly over-determine the system 
and provide a strong verification of the result. In 
addition, we have compared our result with that 
recently reported in Ref. [|| using a completely 
independent method, and find exact agreement. 

4. Hadronic Results 

In Ref. j|, the partonic cross section was 
computed using a series expansion out to order 
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The difference between using that se- 



ries expansion and the exact calculation is quite 
small (less than ~ 1%) as one would expect. In 
light of the intrinsic uncertainty due to scale de- 
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pendence of order ±10%, this difference is com- 
pletely negligible. The convergence of the series 
for scalar Higgs production is shown in Fig. |} 
Preliminary results for pseudo-scalar Higgs pro- 
duction is shown in Fig. || 
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Figure 2. Cross section for pseudo-scalar Higgs 
production at leading order (LO), NLO and 
NNLO (preliminary) j|. The model dependent 
couling g t has been set to unity. 



Figure 1. NNLO cross section for Higgs produc- 
tion using the exact partonic cross section and 
series expansions truncated at the indicated val- 



We obtain our hadronic results by folding the 
partonic cross section with appropriate parton 
distribution functions (PDFs). We use approx- 
imate NNLO parton distribution functions [[[o| 
based on an approximation of the evolution equa- 

tion raa. 
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